Structural and dynamical aspects of avoided-crossing resonances in a 3-level A system 
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In a recent publication [Phys. Rev. A 79, 065602 (2009)] it was shown that an avoided-crossing 
resonance can be defined in different ways, according to level-structural or dynamical aspects, which 
do not coincide in general. Here a simpie 3-level system in a A configuration is discussed, where the 
difference between both definitions of the resonance may be observed. We also discuss the details of 
a proposed experiment to observe this difference, using microwave fields coupling hyperfine magnetic 
sublevels in alkali atoms. 



I. INTRODUCTION 



The concept of "resonance" is ubiquitous in physics but 
it is not always sharply or uniquely defined. The defini- 
tion of a resonance in quantum scattering systems, for 
example, has been subject of endless debates among sup- 
porters of complex plane poles, of a phase shift jump, or 
of other criteria. Typically a resonance implies variations 
of different characteristic properties with respect to one 
parameter within the resonance width, but the extremal 
points for the different properties do not necessarily coin- 
cide; it may even be the case that the variation of a prop- 
erty is completely missing, and the extrema for several 
of them may be shifted with respect to each other, which 
could lead to noticeable errors if the criterion chosen for 
selecting the parameter value is not the most appropriate. 
Resonances appear in summary as multifaccted phenom- 
ena, and a full characterization of their various aspects 
is important to control and optimize specific effects. 

In a recent publication [l| , we have studied the defini- 
tion of a resonance in quantum systems with discrete en- 
ergy levels, in particular those resonances associated with 
avoided crossings. The crossing/avoided crossing sce- 
nario is quite common in many fields of nuclear, atomic, 
or molecular physics such as laser driven trapped ions 
@, El i two level atoms coupled to a cavity mode [3, H| , 
or diamagnetic hydrogen in magnetic fields @. In the 
avoided crossing regions, two eigenvalues of the system 
approach as a parameter of the system A is varied but 
then veer from each other. For a zeroth order Hamilto- 
nian defining the bare levels, the levels do cross at a refer- 
ence value Ao, but a perturbation connecting them causes 
the splitting. The eigenvalues also interchange their char- 
acter so that, in a continuous adiabatic passage following 
one of the eigenvalues through the region, the system 
suffers a significant transformation, being dominated by 
different bare levels on both sides of the crossing. The 
resonance is also characterized by maximal oscillations 
for transition probabilities among the bare levels. As it 
was shown in the parameter values of minimal split- 
ting and of maximal transition probability do not coin- 




FIG. 1. Simple Raman 3-level setup with energy levels 1, 
2 and 3, detunings of the two lasers with respect to atomic 
transitions and coupling strengths (Rabi frequencies) Sli and 

0,2. 



cide in general, defining in this way different, structural 
and dynamical aspects of the resonance. 

In this paper, we propose a simple physical setting, a 
3-level system subjected to a 2-photon transition, where 
this phenomenon may be observed. In Sec. lU the model 
and the different definitions of the resonances are pre- 
sented in detail. The experimental realization is dis- 
cussed in Sec. IIII| and the paper ends with a summary 
and a technical appendix. 



II. THE MODEL 



Consider a 3-level system in a A-configuration (Ra- 
man 2-photon setup, Fig. Q]) which is described, in a 
laser adapted interaction picture, by the time indepen- 
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FIG. 2. (color online) Bare (black solid line) and dressed (red 
dashed line, fii = 0,2 = 0.5^2) energy levels as a function of 
Si/ 6a. 



dent Hamiltonian (h = 1) 



fr = -<y 1 |2)(2| + (<y 2 -* 1 )|3)(3| 

+ %(|2){1| + |1)(2|) 
+ ^(|3)(2| + |2)(3|) 




(1) 



where Oi and O2 are the coupling strengths (Rabi fre- 
quencies) of the different transitions and 5\ and S2 the 
frequency detunings as shown in Fig. [1] When the lasers 
are turned off (Qi = f2 2 =0), the atomic states are un- 
coupled and the energy levels of H cross each other at 
Si = and Si =82- When the coupling lasers are turned 
on, these crossings become avoided crossings and transi- 
tions between the involved atomic energy levels at each 
resonance may occur, see Fig. [2] 

An analytical diagonalization of the full Hamiltonian 
([1]) is possible but the resulting mammoth expressions 
are hardly illuminating. In order to have simple formu- 
lae and gain some understanding about the different as- 
pects of the resonances, approximations will be useful. 
Among the two resonances (avoided crossings) observed 
in the energy spectrum of the 3-level system we shall fo- 
cus on the one at Si = $2- The distance between both 
resonances is S2, and since the energy splitting of each 
avoided crossing is proportional to the Rabi frequencies 
of the coupling lasers fii and f2 2 , the avoided crossings 
will be well isolated (leading to clean transitions) as long 
as Si ~ S2 3> fii, ^2- Under this condition, the state |2) 
is scarcely populated and can be adiabatically eliminated 
to give an effective 2-levcl Hamiltonian as shown below. 



Adiabatic elimination of state |2) and effective 
Hamiltonian 



Substituting the general state \ip) = J2n=i c n(t)\n) 
into the Schrbdingcr equation, the equations for the time 
dependent c n (t) amplitudes are 



ICi 



-C2, 



..fii fi 2 
IC2 = —ci - O1C2 + ~2 C 3, 



IC 3 



fi 2 



C2 + (82 - Si)c 3 . 



(2) 
(3) 
(4) 



The usual, adiabatic-elimination argumenlQ is that when 
Si ~ S2 3> fii,fi2, the population in level 2 remains 
small, nearly zero, and thus 62 (t) ~ 0. In this way, we 
may write 02(f) as a function of Ci and c 3 in the second 
equation and substitute in the other two, 
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(5) 
(6) 



This system corresponds to an effective 2-level Hamilto- 
nian 



Si c 3 . 



He- 




4<5i 

S2 - Si + H 



(7) 



45i wz Wl 1 4.5! 

or, by shifting the zero of energy to make it symmetrical, 



— S e fi fieff 
fieff S e g 



(8) 



which corresponds to an effective coupling of a laser and a 
two-level system with an effective coupling strength f2 ff 
and an effective detuning 5 e s given by 



fieff 



^1^2 

4Si ' 



SeS = r ($2 - Si) 



(9) 
(10) 



As described in when both the diagonal and non- 
diagonal terms in a two-dimensional Hamiltonian depend 
on the same parameter (Si in this simple case), the lo- 
cation of the resonance is not uniquely defined and it is 
possible to use structural and dynamical criteria to define 
the resonance. 



A more accurate and systematic theory, where the exact energy 
levels are obtained by iteration, may be used, see 3 and Ap- 
pendix [X] 
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FIG. 3. (color online) Exact dynamical shift computed by 
numerically diagonalizing the full 3-level Hamiltonian |T|) 
(solid lines) as a function of the ratio between the coupling 
strengths. The approximate expression for the dynamical 
shift (|19[) is plotted with dotted lines. 



B. Structural definition of the resonance 

From the "structural" perspective of the energy-level 
diagram, the resonance may be denned as the point where 
the distance between the two branches of the avoided 
crossing is a minimum. The eigenenergies of H c g are 
easily calculated, 



4<5i 



(11) 



so the minimum distance is given by the condition 



= 0. 




(12) 



This corresponds to a 4 th order equation, whose solution 
will give us the resonance position according to the struc- 
tural criterion. For the condition Si ~ S 2 3> Sli , SI2 , this 
is approximately given by (up to 4 th order terms in the 
frequencies) 



S 2 + 



A8 2 



(f2| 



165| 



46% 



(13) 



C. Dynamical definition of the resonance 

From a dynamical perspective the resonance is defined 
by the value of 8± for which the transition probability 
from state |1) to state |3) is maximum. Using H c g this 
probability is easily computed, 



FIG. 4. Schematic setup with a weak probe field coupling 
states |1) and |3). 



and shows a maximum at <5 e ff = 0, which corresponds to 

(15) 



(Si)d = - [fc + V^ + nf-n? 



n 2 - n 2 (n 2 2 - n 2 ) 



4S 2 



mi 



(16) 



This is also the middle point where the character of each 
dressed energy level changes, as discussed in [l|. The 
cigenstatcs of H e g can be written as 



= sin -|1) + cos -| 3), 



cos-|l) -sin-|3), 



(17) 
(18) 



with tan£> = — f2 e ff/<5eff and real Rabi frequencies fii, 
O2 ■ The change of character of each dressed state is 
centered at the point where the linear combination has 
equal weights (6 = tt/2) for the states |1) and |3). This 
occurs for an effective detuning 5 e g = 0; i.e., this criterion 
coincides with the dynamical definition of the resonance. 

Instead, the expressions (| X3[) and ([To]) do not coincide, 
and are separated by a dynamical shift given in this 
approximation by 



A 



D 



~46f- 



(19) 



which it is plotted in Fig. [3] as a function of the ratio be- 
tween the Rabi frequencies. We see a good coincidence 
between the exact dynamical shift and the approximation 
(|19[) for weak couplings. For strong couplings the pertur- 
bative approach breaks down, and the approximate ex- 
pression deviates from the exact result. In any case, the 
simple form (fTI)|) still gives a good estimate of the effect. 



III. EXPERIMENTAL DETERMINATION 



Pi 



13 



"eff 



s 2 + n 2 



2 

eff 



n 2 



(14) 



The "dynamical resonance" in Eq. (|16p can be eas- 
ily determined experimentally by preparing the system 
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in state |1) for each <5i and looking for the maximum 
probability of finding state |3). The experimental de- 
termination of the minimum level splitting ("structural 
resonance" ) requires some more work. One way is to use 
a third, auxiliary, weak probe field connecting states |1) 
and |3) as depicted in Fig. [4j This transition may be an 
electric-dipole forbidden transition, such as a magnctic- 
dipole-allowed transition. Magnctic-dipole transitions 
are usually much weaker than electric-dipole transitions, 
a good thing in this context since we are interested in 
probing the dressed energy levels without excessively per- 
turbing the original system. 

The Hamiltonian describing the full system (including 
the probe field) takes the time-dependent form 



H(t)=H + W(t), 



(20) 



since, in general, there is no interaction picture in which 
the full Hamiltonian is time independent. Here, H is the 
Hamiltonian of the original system already given in Eq. 
([!]), and the time-dependent perturbation is given by 



W(t) = ^ (|3)<l|e 4 
v = Si — 62 — S 



H.c) 



13 j 



(21) 
(22) 



<5i3 being the detuning of the probe field with the 1) o 
3) transition, see Fig. SJ 

We shall examine hereafter the resonance at Si ~ S2, 
see Fig. [21 The dressed states will be labeled with in- 
creasing energy (ei < €2 < £3), so we have to measure 
the energy difference between e 3 and £2 for determining 
the structural resonance. 

We shall now consider if as a zeroth order Hamiltonian 
weakly perturbed by W(t), and use (time-dependent) 
perturbation theory to obtain the transition rate from 
dressed state \e2) to dressed state (£3), 
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Ae 2 



(cos 2 ut — cos vt cos Act) 



, (23) 



with a-, 



£3 ~ £2- P\e 2 )^\e 3 ) 



*ij - <e 3 1«) 0' |e 2 > and Ae(£i) 
will show peaks at v « ±Ae. Thus, by changing the 
probe detuning Siz (sweeping the value of v) and mea- 
suring the corresponding transition rate for a fixed set of 
parameters of the probeless system, the energy splitting 
between levels £2 and £3 is determined. Following the 
same procedure for different values of cq, it is possible 
to find the minimum splitting and identify the structural 
resonance. 

Note that the positions of the maxima of P| e2 )_>.| e3 } 
will only be located exactly at v = ±Ae in the long 
time where when Sln ((^^)*/ 2 ) ^(z, ± Ae) and the 
contribution of the crossed term in Eq. (|23[) becomes 




[(5 1 ) s exact -(8 1 ) s probe ]/5 ; 

V 8 2 
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FIG. 5. (color online), (a) Exact energy splitting Ae = £3 — 62 
obtained by diagonalizing the full 3-level Hamiltonian given 
in Eq. (fl| (black-solid line), compared to the energy splitting 
obtained by identifying one of the maxima (in this case, the 
maxima around v ~ Ae) of the P| E2 )_>| C3 ) transition proba- 
bility (red-dashed). In the t — > 00 limit both lines converge 
(the calculation is done for S-iljl-E = 125). (b) Difference be- 
tween the exact (probeless) position of the minimum splitting 
(structural resonance) and the structural resonance obtained 
from the probed system as a function of time (black-solid 
line). This difference goes to zero in the long time limit as 
expected. The red-dashed line corresponds to the exact value 
of the dynamical shift A o (which is indeed the precession re- 
quired in the measurement of the structural resonance) . The 
time marked by an x in (b) corresponds to t — 2tt/Ad, the 
lower bound of the time requirement in order to resolve the 
dynamical shift, see Eq. ([24)) . Sli = 0.2S 2 , 0,2 = 0.552. 



negligible. At short times, the positions of the maxima 
arc shifted due to the <5i dependence of the 

Gt<^2 S ^ SGG 

Fig. EK- To resolve the dynamical shift by this method, 
this effect should be smaller than the dynamical shift 
itself, which is indeed achieved at sufficiently long times, 
as shown in Fig. [SJ). 



As time increases the peaks of Pi 



become nar- 



rower (the width of each peak goes like 2ir/t), so in order 
to be able to resolve the dynamical shift, the probe beam 
should be applied for a time satisfying 



t > 



2tt 

a7 



(24) 



In summary, for large enough times both effects (the shift 
due to the <5i dependence of the o^ 's and the width of 
the peaks in order to resolve the dynamical shift) can be 
overcome. Note also that condition (|24[) . which ensures 
narrow peaks, is more demanding than the times required 
to get rid of the shift due to the Si dependence of the 
Oij's, see Fig. [5b. 

The height of the peaks grows with time as ~ £l 2 t 2 /4, 
so to keep the perturbative treatment valid, this maxi- 
mum probability has to be smaller than one (weak probe 
field). Combining this low probe intensity condition with 
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the long-time condition given above, we end up with a 
condition for the probe- field amplitude f2 p , 



(25) 



Actually this is just an upper bound. The exact growth 
of the height with time is given by (a^,) 2 £2 2 i 2 /4, but the 

values of the matrix elements a^ k arc bounded between 
and 1. 



A. Discussion of specific systems 

The most obvious setting for a Raman-transition ex- 
periment is driving optical stimulated Raman transitions 
in alkali atoms. Unfortunately, this appears to be a dif- 
ficult scenario in which to study this effect. Taking 87 Rb 
as an example, for driving stimulated Raman transitions 
between hyperfine ground levels, using lasers nearly res- 
onant with the D2 line (5 2 Si/2 — > 5 2 P3/2 transition), 
typical parameters are a detuning 5\ as 82 = 2ir ■ 10 GHz 
and Rabi frequencies f2i = SI2 = 2ir ■ 200 MHz. These 
parameters give a lowest-order dynamical shift [Eq. ([TTj)) ] 
of 400 Hz. However, with an excited-state decay rate 
of r = 2ir ■ 6.1 MHz, the rate of spontaneous scattering 
from the Raman fields is around R sc « r(f2 2 + ri|)/8(5 2 , 
or about 3.8 kHz. The problem here is that the dressed 
states will be broadened at the kHz level, and the in- 
teraction time of the probe will be limited, so that the 
the resolution of the probe will be too poor to resolve 
the dynamical shift. Decreasing the scattering rate also 
does not help much; for example, increasing the detun- 
ing to 100 GHz leads to a scattering rate of only 38 Hz, 
but a dynamical shift of only 400 mHz. The scattering 
rate becomes comparable to the dynamical shift for a de- 
tuning of only 1 GHz, which is realistically too small for 
precision measurements. 

A more promising experimental realization is pos- 
sible by driving microwave transitions in the hyper- 
fine structure of the ground electronic level of atoms. 
Here, spontaneous emission is completely ignorable, as 
the magnetic-dipolc transition lifetimes are much longer 
than any reasonable laboratory time scale. In par- 
ticular, we consider here the n 2 Si/2 ground state of 
alkali atoms, which is split into two hyperfine levels, 
F = I ± 1/2, where / is the nuclear-spin quantum 
number. The three hyperfine sublevels corresponding to 
the setup in Fig. Q] arc |1) = \F = I — 1/2, m F = —1) 
and |3) = \F = I + 1/2, m F = —1) for the two Raman- 
coupled states, and |2) = \F = I + 1/2, m F = 0) for the 
intermediate ("excited") state. The degeneracy of the |2) 
and 1 3) states is broken by applying a magnetic bias field 



-£*bias — 



X AE, 



his 



Mb (5j - .9/)' 



(26) 



where A£hf s is the zero-field hyperfine splitting, g, and 
g 1 are the electronic and nuclear (/-factors, respectively, 



and ^ := (J+ 1/2) -1 . This represents the center of an 
avoided crossing of the |1) and |3) states, and thus the 
splitting at this bias-field strength, 



AE; 



31 



VT 



X' 



AE 



ill's- 



(27) 



is insensitive to first order to bias-field fluctuations. This 
reduces the need for stringent experimental control over 
magnetic fields, and reduces the most important system- 
atic error in measuring the Raman resonances. With the 
same magnetic field, the energy of the "excited" |2) state 
is above that of the |3) state by an amount 



Ai?23(-E>bias) = S/^B-Bbi, 
1 AiiMs r r— 



XT 



(28) 



For example, for 87 Rb, with AE Ms = h ■ 6.835 GHz and 
I = 3/2, the bias field is -Bbias = 1-219 kG, and the split- 
tings are AE31 = h- 5.919 GHz for the (nominal) Raman 
resonance, and AE23 — h ■ 860 MHz for the ft.2 driv- 
ing transition. The remaining (fii) driving transition is 
given by the sum of the other two transition frequencies, 
or AE21 = h ■ 6.779 GHz. Both Raman driving transi- 
tions are driven by circularly polarized fields, while the 
probe field is driven by linearly polarized field. 

Continuing with the 87 Rb example, the Raman fields 
may be applied with Rabi frequencies of Qi = £7 2 = 
2ir ■ 300 kHz, corresponding to field intensities of about 
7.6 W/cm 2 on both transitions. Microwave fields of this 
intensity, for example, have been realized around 6.8 GHz 
in the near-field of an atom chip to manipulate a Bosc- 
Einstein condensate of 87 Rb [8[. Thus, a field of this 
strength for the 6.8 GHz transition is feasible, and the 
field for the 860 MHz transition should similarly pose no 
problem. For a Raman detuning 5± ~ 8 2 — 2ix ■ 1 MHz, 
the lowest-order dynamical shift from Eq. (|19l) is 2.0 kHz. 

In the choice of parameters here, it is also convenient 
to have very different Raman transition frequencies (6.8 
and 0.9 GHz for the fii and fl2 fields, respectively), to 
control the secondary ac Stark shifts that we have not ex- 
plicitly accounted for. That is, for example, the 6.8 GHz 
Oi field driving the |1) — > |2) transition also couples 
the the |3) — > |2) transition at 0.9 GHz, albeit much 
farther off resonance. As long as 62 is held fixed, the 
Stark shift of |1) due to the SI2 field is inconsequential, 
as it simply causes a common shift of both structural and 
dynamic resonances. However, the Stark shift of \2) due 
to the fii field depends on 5\, and thus can cause an addi- 
tional contribution to the dynamical shift Aj> However, 
this effect is suppressed by the ratio of the detuning Si 
from the |1) — > |2) transition to the detuning from the 
|3) — > |2) transition. This effect should thus be smaller 
than the lowest-order shift of 2.0 kHz by a factor of about 
10~ 4 , and is therefore negligible. Note also that it is im- 
portant to have Raman detunings much smaller than the 
transition frequencies, in order to suppress the effects of 
Bloch-Siegert shifts. By a similar argument, the contri- 
bution of the Bloch Sicgcrt shifts should be of the same 
order as the secondary ac Stark shifts. 
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Note that uncertainties in the microwave frequencies 
are negligible on the scale of kHz, so long as the fields are 
derived from digital synthesizers. However, the splitting 
at each detuning must be determined to an accuracy finer 
than the 2.0 kHz shift. Thus, to resolve this shift of 
2.0 kHz, the (6.8 GHz) probe beam should be applied 
for a time much longer than 500 ixs [Eq. (|2"4")) ]. with a 
Rabi frequency small compared to 2ir • 1 kHz [Eq. (|25p]. 
The probe field then requires a correspondingly much 
lower intensity, as compared to the Raman fields. The 
atoms will also need to be well-confined on ms time scales, 
without inducing spontaneous emission. Loading laser- 
cooled atoms into a dipole trap — formed by the focused 
light of a CO2 laser — accomplishes this confinement, with 
negligible perturbation to the hyperfine structure of the 
ground electronic state. 

Care must also be taken in preparing the atoms for 
the probe measurement. Since the goal is to measure the 
splitting of the dressed states at a particular detuning, as 
described above, we must prepare the atoms in only one 
of the dressed states. This is effected, for example, by 
first optically pumping the atoms (in the absence of the 
Raman fields, and with only a small magnetic bias field 
of the order of 100 mG to prevent mixing of states) into 
the |1) = \F = l,m F = —1) bare state. This is accom- 
plished by driving the 5 2 S 1/2 , F = 1 — >5 2 P 3/2 , F 1 = 1 



optical transition with circularly polarized light, while 
optically depumping the F — 2 ground hyperfine level. 
The 1.2 kG field should then be turned on adiabatically 
to produce the correct level configuration without induc- 
ing any transitions. The Raman fields should then also be 
turned on adiabatically, but far from Raman resonance. 
They can then be adiabatically chirped to the desired 
Raman detuning, transferring the atoms from |1) to 1 62) - 
The probe field should then be activated to attempt to 
drive atoms to the other dressed state |e3). Finally, the 
Raman fields should again be detuned and adiabatically 
turned off, and the magnetic field turned off adiabatically 
as well. The population transferred to ^3), and thus to 
3) = \F = 2,m F = —1), is then measured by fluores- 
cence detection of the F — 2 population. Finer resolu- 
tion of the Raman splitting is also possible by employing 
a Ramsey-interference technique, applying the probe in 
two pulses separated in time. 



IV. SUMMARY 

As discussed in avoided crossing resonances be- 
tween dressed energy levels are not uniquely defined. In 
this paper a 3-level atom in a A configuration has been 
proposed as a simple physical system where the distinc- 
tion between structural and dynamical aspects of a res- 
onance could be observed. By adiabatically eliminating 
the third level, an effective 2-dimensional Hamiltonian 
has been obtained, from which an approximate analyti- 
cal expression for the dynamical shift (between structural 
and dynamical resonances) has been given. 



While the dynamical resonance is in principle easy 
to observe by measuring the maximum rate of a given 
atomic transition, the determination of the structural 
resonance is more delicate. We have proposed a method 
consisting of a weak probe field. In order to resolve the 
dynamical shift by this method, a low-intensity probe 
field must be applied for sufficiently long times. An ideal 
setting is provided by microwave transitions in the hyper- 
fine structure of the ground state level of alkali atoms. 
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Appendix A: Resolvent method 

It is possible to improve systematically the adiabatic 
approximation by using a more accurate method, where 
the energy levels are given exactly by an implicit Hamil- 
tonian, [J, [BJ . In order to use this resolvent method, we 
may divide the starting Hamiltonian ([1]) as H = Hq + V 
with 




H = I -5j 

6 2 - <Si 



1 



fti 



v = - ^1 n 2 
n 2 



(Al) 
(A2) 



The bare energy levels (eigenergies of Hq) corresponding 
to the bare states |1), |2) and |3) are given respectively 
by 



-1 
,(0) 



0. 



Si, 



^=82-6,. 



(A3) 
(A4) 
(A5) 



Around Si « 5 2 , energy levels and £3 are degenerate 
but will form an avoided crossing when the lasers are 
turned on. Around this value of 8\, the system will be 
described by an implicit effective 2D Hamiltonian Q, 



H. 



eff 



-^31 



Rll i?i 3 

£ (0) 
e 3 



R 



(A6) 



33 



where Rij — (i\R\j) are the matrix elements of the level 
shift operator R, 

00 

R{E) = Y. PV 



71=0 



Q 



E-Hn 



V) P, 



(A7) 
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with P = |1}(1| + |3> <3 1 and Q = 1 - P = |2)(2|. The ef- 
fective Hamiltonian (|A6|) can be written in a symmetrical 
way by changing the zero of the energy 



with 



Sett 

c 



—S e e R13 

R31 S e ff 



c 



1 
1 



(A8) 

(A9) 
(A10) 



Explicit expressions of the needed elements are exactly 
given by 



R11 

R33 
R13 



4(£ + <5i) : 
4(£ + c5i)' 
4(S + (5i) 



#31 1 



(All) 
(A12) 
(A13) 



where we have taken into account that e 



(0) 



-Si, and 



C 



Si 



Si 



n? 



4(£? + <Ji) 



4(£? + *i) 



(A14) 
(A15) 



Note that the elements of the implicit Hamiltonian de- 
pend on the eigenenergy E. Thus the energy levels will be 
obtained by iteration. Once the energy levels are calcu- 
lated (with a given precision) the different resonance loci 
will be calculated as described above. Choosing E = 
for the first iteration corresponds exactly to the adiabatic 
elimination approximation in Sec. |TTJ 
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